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We develop a hydrodynamic description of the collective modes of interacting liquids in a quasi- 
one-dimensional confining potential. By solving Navier-Stokes equations we determine analytically 
excitation spectrum of sloshing oscillations. For parabolic confinement, the lowest frequency eigen¬ 
mode is not renormalized by interactions and is protected from decay by the Kohn’s theorem, which 
states that center of mass motion decouples from internal dynamics. We find that the combined 
effect of potential anharmonicity and interactions results in the depolarization shift and final lifetime 
of the Kohn mode. All other excited modes of sloshing oscillations thermalize with the parametri¬ 
cally faster rate. Our results are significant for the interpretation of recent experiments with trapped 
Fermi gases observed weak violation of the Kohn theorem. 

PACS numbers: 67.10.Jn, 72.15.Lh, 72.15.Nj 


Introduction. Properties of quantum liquids in one- 
dimension (ID), as realized experimentally in nanoscale 
semiconducting wires, carbon nanotubes, laser traps of 
cold atoms, and edge channels of the quantum Hall effect, 
continue to attract tremendous attention in the physics 
research (see Refs.^i®! for recent reviews and references 
herein). With increasing sophistication in high precision 
measurements and techniques these systems provide se¬ 
rious tests for the existing theoretical models, such as 
Luttinger liquid theorj®^, and ultimately challenge their 
completeness. For example, the powerful approach of the 
Luttinger liquid formalism allows to account for the in¬ 
teraction effects nonperturbatively. However, this model 
does not adequately describe relaxation phenomena due 
to built-in approximation of the linearized quasiparticle 
dispersion, which by virtue of the kinematic constraints 
effectively closes the phase space available for inelastic 
scattering. In certain special cases, the lack of relaxation 
may be a generic property of the many-body system be¬ 
cause of its complete integrabilitj®^. Alternatively, van¬ 
ishing relaxation rates may happ en be cause of the reasons 
prescribed by the Kohn theorenj- ^° * ^^ l 

Motivated by recent experimentJ^^HUl study re¬ 
laxation of collective excitations in interacting two- 
dimensional (2D) systems confined along one of the 
two spacial dimensions. These systems interpolate be¬ 
tween strictly ID limit of Luttinger liquids and the two- 
dimensional Fermi liquids. The geometrical confinement 
in such systems splits the single band spectrum into mul¬ 
tiple ID subbands. The convenient and practically jus- 
tihed model idealization is the interacting particles con¬ 
fined by a harmonic potential. In this case the ID sub¬ 
bands are equidistantly separated by a frequency of os¬ 
cillations uj± across the channel. Similar to the spectrum 
linearization in the strictly ID liquids, harmonic approx¬ 
imation in the quasi-ID case on one hand simplifies the 
dynamics, and at the same time does not allow for proper 
description of thermalization processes. One necessarily 
has to account for the confinement anharmonicity, which 
thermalizes the motion across the channel in much the 


same way as the spectrum nonlinearity causes the re¬ 
laxation of charge and spin excitations in ID quantum 

^irepmSI 

Despite the similarity with ID, the relaxation of 
transversal excitations has a few distinct features setting 
these two problems apart. The kinematical constrains of 
momentum and energy conservations operational in ID 
are less restrictive in quasi-lD. In contrast to the ID 
case, which require three-particle scattering processes, 
the two-body collisions do cause the relaxation via the 
inter-subband transitions. Thermalization in quasi-lD 
may nevertheless be prohibited due to the Kohn theo¬ 
rem rather than kinematical restrictions. 

This theorem states that the motion of the system as 
a whole is unaffected by interactions. Classically, it fol¬ 
lows as the translationally invariant interaction energy is 
insensitive to the system displacement as a whole. For 
the same reason, quantum mechanically, the interaction 
drops out o f the center of mass Heisenberg equation of 
motioilldJ211 jn a quantum Fermi liquid the Kohn theo¬ 
rem follows from the Landau Galilean invariance relation 
between the quasiparticle effective mass and the first an¬ 
gular harmonic of the interaction amplitud^^. 

In all of the above cases the Kohn theorem states that if 
the confining potential is harmonic the collective sloshing 
oscillations proceed without decay. The frequency of the 
Kohn, or so called sloshing mode, uj±, is ins ensitive to 
interaction, temperature and particle statisticJ^^^^IIall 

This fact makes the observation of the Kohn mode pos¬ 
sible in a wide variety of systems. In semiconductor quan¬ 
tum wires the Kohn mode is observed in optical trans¬ 
mission at far infrarecP^lMI, jn trapped ultracold Fermi 
gas of ®Li the Kohn mode of a half KHz frequency was 
excited by sudden displacement of the tr ap and detected 
by absorption imaging of a released cloucP^ESl, 

The weak violation of the Kohn theorem due to an¬ 
harmonicity observed in the above three classes of sys¬ 
tems plays a key role in the data interpretation. In 
the case of the semiconductor quantum wires it controls 
the line broadening and the higher harmonics of light 
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FIG. 1. (color online) The definition of the displacement 
field (j>{z, t) in Lagrangian formnlation. At the time, shown 
as vertical axis, t = 0 the flnid is contained in the interval 
\z\ < a. As time progresses the particles of a flnid move as 
indicated by narrow arrows pointing upward. The particle 
located at 2 : at t = 0 is shifted to the new position 2 + (j){z, t). 
A fluid volume shown as a narrow (red) horizontal rectan¬ 
gle occupies a segment [ 2 , 2 -I- dz\ at time t = 0. At a later 
time t > 0 this volume is displaced and occupies the segment 
[2 -I- 4>{z, t), z + dz + (j>[z + dz, t)]. The fluid volume expands 
(shrinks) if d4>{z,t)/dz > (<)0. The expansion (shrinkage) 
translates in the decrease (increase) in the density respec¬ 
tively. 


transmissiorP^. The observed sloshing frequency of an 
atomic cloud in the optical trap shows systematic devi¬ 
ations from the Kohn theorem predictiorP^. Such devi¬ 
ations grow with heating as the expanding atomic cloud 
senses a progressively less parabolic confining potential. 
Here we concentrate on two fundamental aspects of Kohn 
theorem violation: (i) depolarization shift of the slosh¬ 
ing frequency, and (ii) final lifetime of sloshing oscilla¬ 
tions. We approach this problem based on a very general 
grounds of hydrodynamic theory, which accurately de¬ 
scribes most liquids at length scales long compared to 
the particle-particle mean-free path. 

Hydrodynamic theory. A hydrodynamic description is 
based on the existence of slow variables associated with 
locally conserved quantities such as number of particles, 
momentum and energy. The motion of the liquid is de¬ 
scribed by the Navier-Stokes equations which in Eulerian 
continuous field coordinates can be put in the fornP^ 

dt{pvj) =-diliij - pdjU, ( 1 ) 

that guaranties the momentum conservation, which holds 
in ideal and nonideal liquids alike in the absence of con¬ 
fining potential U. The stress tensor of a two-dimensional 
fluid giving rise to the Navier-Stokes equation iP^ 

Hy = S^jP + pviVj - (SijdkVk - T]{diVj + djVi - SijdkVk). 

( 2 ) 

Here rj, ( are the first (shear) and second (bulk) viscosi¬ 
ties, and P is pressure. In the following we focus on a 
ID flow in 2 -direction of a two-dimensional liquid occu¬ 
pying the strip \z\ < a, so that velocity vector field can 
be taken in the form v = v{z,t)ez and Eq. ([^ simplifies 


to 


p{dtv -I- vdzv) = -dzP - pdzU + dz{ridzv). (3) 

When writing this equation we made use of the con¬ 
tinuity equation dtp + dz{pv) = 0 and employed stan¬ 
dard assumption ij ^ For the purpose of our study 
it will be convenient to use a particle description of 
Navier-Stokes equation. In this approach the coordi¬ 
nate 2 labels an equilibrium position of a fluid parti¬ 
cle, while its location at later time t, z + 4>{z,t) defines 
the displacement field (j){z,t), see Eig. The density 
is p{z + 4>{z,t),t) = pt)(z)/{l P dz(j){z,ty), the velocity 
v{z + (j){z,t)) = The linearization of Eq. @ 

is equivalent to the linearization in (j). To the leading 
order, 

p = po + Sp, Sp=-{pocj)y , v = (j), (4) 


where po is a stationary equilibrium density distribution, 
the notations /' = df jdz and / = df /dt are introduced 
and the pair of arguments ( 2 , t) common to all the func¬ 
tions is omitted. The parametrization (|^ of Jp and v by 
a single displacement field automatically satisfies the lin¬ 
earized continuity equation, 5p-|-(po0)' = 0. The concept 
of the displacement field is further illustrated in Fig. 
For the solutions of the form (j){z, t) = the equa¬ 

tion ^ in the parameterization, Q reads, see App. 

^'^X =-vlx!' + U'x' -iuj [ —(zzpox")': (5) 

J Zq PO 

where y{z) = rj{z)/p{){z) is the kinematic viscosity, and 
Vs = y dP/ dpo (6) 


has a meaning of a local speed of sound that depends on 
2 only through the equilibrium density po. We show that 
the results are independent on arbitrary zq. 

The spectrum of collective excitations and their decay 
rates as given by Eq. ([^ depend on the details of the 
confining potential and dependence of the viscosity on 
density. Here for definiteness we consider the confining 
potential per unit mass with a weak quartic anharmonic- 

ity, 


W I 2 62’ 


( 7 ) 


where m is the mass of an individual particles, and the 
constant SU cx e is added in such a way that the spatial 
extent of the liquid stays the same as for e = 0. In other 
words 6U -fea'^/dm = 0. This choice is not obligatory, yet 
convenient in the further calculations. The anharmonic- 
ity oc 6 modifies both the first and the second term in the 
right-hand-side of Eq. Indeed, ^^( 2 ) = 2TTp{z)jm^ = 
w^(a^ — z'^)j2, valid for the parabolic confinement, ac¬ 
quires a correction 5v1{z) = —e{z'^ — a^/ Am. With these 
observations in mind we multiply Eq. p|) by 2/, rescale 
coordinates z ^ zja and introduce X^= 2a;^/w^ to find 

A^X + (1 - ^')x" - 2zx' + Vs + V, = 0. (8) 
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The perturbation term due to anharmonicity can be cast 
in the form 


= - 


ea 


2mu)\ 




(9) 


The generalization of Eq. ([^ to arbitrary confining po¬ 
tential is accomplished by replacing [{z"^ — l)x']' —>■ 
[{U{z) — U{a))x']'■ Consequently, our results for the de¬ 
polarization shift are straightforward to modify for arbi¬ 
trary confinement. We emphasize that K, is Hermitian. 
This property guaranties that the anharmonicity alone 
causes only the frequency shift but no dissipation. We ar¬ 
gue below that only the combination of the anharmonic¬ 
ity and viscosity leads to the dissipation of the Kohn 
mode. The exact z-dependence of the viscosity v{z) is 
specific to the model. This however only influences the 
numerical prefactors in the final results and we take the 
expression for the viscosity from the theory of Fermi liq¬ 
uids, rj oc VFP^^ where mean fee path is ^ oc ypEp/T'^ 
with vp and Ep being Fermi velocity and energy, respec¬ 
tively. For this case ^{z) = Cp'^{z)/m^T^ where C is the 
numerical factor of the order of unitj®!. In the above 
specified dimensionless notations this results in 

K = -^X^B r [(1 - (10) 

Jzo ^ ^ 


where we have introduced B = a^ujp with vq = 

v{z = 0). We proceed with the analysis of Eq. (|^. 

Results. As the first step let us discuss the eigenmodes 
of an ideal fluid confined to a harmonic trap. For that 
purpose we neglect anharmonicity and interaction effects 
implicit in the viscosity term of Eq. (|^, i.e. we set 
Ve = Vrj = 0. What remains is familiar Legendre equa¬ 
tion and we therefore immediately read off its solutions 
= n{n + 1) with n = 0,1, 2,... so that 
the eigenfrequencies are 


ln{n + l) 

^-• ( 11 ) 

The n = 0 gives an equilibrium since Xn=o = const 
this velocity is identically zero. The n = 1 is a Kohn 
mode wi = u}±. We thus found the whole hierarchy of 
eigenoscillations, they are Legendre polynomials 

Xn(^) = (12) 

and the related velocity fields are u„(z) = d^Pn- Re¬ 
markably, the same spectrum of collective modes as (0 
was obtained recently for the longitudinal oscillations of 
the ID Coulomb chain^^. In this systems the collective 
behavior sets in due to the long range Coulomb forces 
rather than collisions. 

Next we discuss the significance of perturbation terms 
on the spectrum of collective modes. Let us first con¬ 
sider anharmonicity alone, i.e. we set = 0 and 14 ^ 0 


in Eq. (|^. As 14 is Hermitian the spectrum remains 
real. As a result all the eigenmodes remain undamped as 
expected in the absence of collisions. To find the depo¬ 
larization shift of the Kohn mode, to the leading 

order in e we write = 2 + where the cor¬ 

rection term is found from the first order perturbation 
theory with unperturbed solutions given by Eq. (12), 




The matrix element gives a factor of 8/5 which even¬ 
tually translates into the correction to the Kohn mode 
frequency 




3ea^ 

lOmwj, 


(14) 


We have checked that Eq. (14) agrees with the result ob¬ 
tained by the method of moments suggested in Ref.^^. 
Note however that the depolarization shift in the col- 


lisioneless regime differ from Eq. (14) by a nonuniver- 


sal numerical prefactor. For instance, for the contact 
interaction we obtain by direct perturbation theory in 
Fermions = 3ea^/5ma;_L,l^. Corrections to other 

eigenfrequencies can be computed in the same fashion 
with the result ( |l4| ) up to a numerical coefficient. 

Before considering the generic case, it is instructive 
to verify the Kohn’s theorem within our hydrodynamic 
approach. It amounts to the statement that Xi(z) oc z 
remains the solution of Eq. (|^ with the frequency wi = 
ujp even for 14; 0 provided only 14 = 0. This in turn 

will be proven once we show that for any non-negative 
n, {Pi\Vr,\Pn) = {PnK\Pi) = 0. Clearly KJPi) = 0 as 
[Pi(z)]" = 0. On the other hand from Eq. (|^ it follows 
that 


{Pl\Vu\Pn)o^j dzPi{z)J 


^ dz 
,Po{z) 


Mz)po{z)[PM)]'y- (15) 


Realizing that Pi(z) = —dzPo{z)l2 and integrating by 
parts one concludes that {Pi\V,z\Pn) = 0 for all n thus 
proving the Kohn theorem in the present context. 

It follows that the Kohn mode may acquire a finite life¬ 
time only when the anharmonicity is included. Yet the 
perturbation 14 is Hermitian and by itself is insufficient. 
We therefore consider both perturbations, and look for 
the corrections that are first order in each of the two. In 
the second order perturbation theory such a correction is 
of the form 


5(2)a 


2 

CJl 


oo 

-iXuiB ^ 
n—1 


y^V(2n + l)/2 

2 — n{n + 1) 


[{Pl\V,\Pn){PnmPl) + {PimPn){PnK\Pl)]. (16) 


As we saw above both terms in this equation are zero as 
the Kohn mode does not couple to any other mode by a 
viscosity term 14,, and = 0. 
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Inevitably we have to consider the third order correc¬ 
tions to Atjj. The third order correction to the energy Ei 
of a state |Z) is 


5(3) = 


imni\ {l\V\m){m\V\l) 

{Err^-EiY • 


(17) 


To apply this expression to our problem we identify 
11) = Pi as a Kohn mode and V = 14-1-14. We ob¬ 
serve that for the perturbation terms given by Eqs. ^ 
and ( 10 ) the following properties of matrix elements hold 
(Pill^n) = (PniKiPi) = {Pi\V,\Pi) = 0 for all n. Fur¬ 
thermore, for the quartic anharmonicity under consider¬ 
ation the only nonzero off-diagonal matrix elements are 
(P 1 II 4 IP 3 ) = (P 3 II 4 IP 1 ) with the rest of matrix elements 
(Pn|14|Pi) = 0 for n ^ 1,3. We thus have, accounting 
for all the normalization factors of eigenoscillation modes 

pi) 


J(3)A2 = _ — 

8(A2^-AS3) 


2(Pl|F.|^3)(P3|K|P3)(P3|f4|Pl). 

(18) 

Using the explicit expressions ([^ and (10) we find for 
the matrix elements 


dzPs{z). 


dz 


’-1 

»i 


j ^dzp,{z){{z^-i)mz)]'}'= 1, 


/-I 

and eventually 


^_^,-{(i-.^)3[P3(z)rr = -^,(i9) 


( 20 ) 


j(3)\2 28 a/2 Up / ea^ 

“ 125 [mujl 


( 21 ) 


This result enables us to hnd imaginary part of the Kohn 
mode wi = u}j_ + 5a;i -I- which corresponds to its 
attenuation with the rate 


2 9 

__i r'oe a 

U — 2 4 


Suji 

UJ± 




( 22 ) 


where we omitted numerical factors of order unity for 
brevity. This expression constitutes the main result of 
our work and has straightforward interpretation. The 
higher excitations modes not protected by the Kohn the¬ 
orem decay with the rate ^ v^jc?. As vq has dimensions 
of the diffusion coefficient, this is the typical rate of the 
momentum relaxation. The ratio {Suji/uj±)'^ is the prob¬ 
ability of finding the system in higher modes. 

Discussions. Hydrodynamic description requires short 
equilibration length £. Thus validity of our theory is 
limited by the condition i a, which imposes certain 
restriction on temperature. Specifically, for the Fermi 
liquids £ = vpTee is determined by collisions with the 
typical rate ^ T'^/Ep. Since a;j_ ~ vp/a, hydrody¬ 
namic regime is realized at temperatures T > Th above 


the crossover scale Th ^ y/ujjTTp ^ Ep/y/N, where N 
is the number of occupied sub-bands of the transver¬ 
sal quantization. It also follows that with necessity hy¬ 
drodynamics requires T ^ ujp. While this inequality 
is reasonably satisfied for the cold gases that are con¬ 
fined by a very shallow potential, it obviously breaks in 
the ultra-cold limit where collisionless regime prevails. 
In the latter case attenuation coefficient of the Kohn 
mode is expected to follow quadratic temperature depen¬ 
dence oc aT'^/Ep based on the Pauli principle and 
phase space restrictions argument, whereas in the hydro¬ 
dynamic regime oc l/T^ in accordance with Eq. (22 1 . 
The nonmonotonic temperature dependence of the decay 
rate has been observed experimentalljES. 

Our hydrodynamic approach has interesting parallels 
with the Luttinger liquid description of collective modes 
in confined inhomogeneous ID gases 1^. The eigenvalue 
equation for the normal eigenmodes in that case, analo¬ 
gous to our Eq. ([^, is given by 


-^lxn{z) = v{z)K{z)d^ 


K{z) 


C>zXn{z) 


(23) 


where Luttinger liquid interaction parameter satishes the 
relation v{z)K{z) = 7 rp(z)/m^. This equation is supple¬ 
mented by the boundary condition y„(±a) = 0 and nor¬ 
malization condition dzxj(z)xj(z)/v(z)K(z) = 5^. 

For the particular choice of v(z) = vo\/l — z^ja^ and 
K{z) = Kq{1 — the solutions Xn(z) are obtained 

in terms of Gegenbauer polynomials with the spectrum 
of excitations = (r;o/a)^(n -|- l)(n -I- 27 -|- ijMiH, jn 
the model of 7 = 2 , the problem simplifies to the case of 
Legendre polynomials^ with the spectrum of excitations 
analogous to our result ( 0 . Another interesting limit 
is 7 = 0, which corresponds to the case of the Tonks- 
Girardeau gas, where the Gegenbauer polynomials re¬ 
duce to Chebyshev polynomials. Inclusion of dissipative 
terms into Eq. (23) requires consideration of corrections 


to Luttinger liquid model which account for the inelastic 
scattering of bosons and ultimately describe equilibra¬ 
tion processes. As recently shown such generalization is 
possible both in the limit of wealP^ and stron^^ inter¬ 
actions and application of this formalism to the problem 
of decay of collective modes is an interesting problem for 
future research. Along this rout one may hope to find a 
unihed description, which interpolates between the quan- 
turrP^ and classicaP^ hydrodynamic regimes of Luttinger 
liquids, and which is broadly applicable for arbitrarily 
strong interactions. 
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Appendix A: Derivation of equation ([^ 


Writing 


Here we derive the Eq. ([^ of the main text. In the 
parametrization (11> the left hand side of Eq. @ takes 
the form, 


p{dtv + vdzv) « • 


(Al) 


To linearize the right hand side of Eq. ([^ we note that 
the pressure is hxed by the density via the equation of 
state such that 

P{z,t) =P[p[z,t)] K P[pq[z)\ - vl[po{z)]{pQ(l))' , (A2) 

where the velocity Vg is defined in Eq. §. At equilib¬ 
rium, Eq. @ yields 

v^p', = -poU' (A3) 

We have therefore, 

-P' - pU' « [vUpocl^yr + {P0<pyu' (A4) 


Writing {p^cfyu' = [p[j(l)U']' — po4>U'' and using (A3) we 
obtain, 


KVo<^T = PoK^</>T + Pob«</>'] 


and using (A3) again we obtain 

-P' - pU' « po - PO^'U' - Pocj^U" . (A6) 

The third, viscosity term on the right hand side of Eq. ([^ 
reads 


dzivd^v) = [rj^'Y 


(A7) 


Substituting Eqs. (All, (A6) and (A7) in Eq. ^ of the 
main text we obtain 

- 4>'U' - cj^U" + pY^ir,^']' (A8) 

For the solutions of the form (j){z, t) = we obtain 

the equation, 

= [vlx!']' - X"U' - x'U" + (-*a;)po'M' 

(A9) 


-P' - pU’ « [VgPoft)']' - {Po<y)U'' (A5) Integration of Eq. ^ over z yields Eq. 
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